Preprint typeset in JHEP style - HYPER VERSION 



Minkowski solution of Dyson-Schwinger equations in 
momentum subtraction scheme 



Vladinn'r Sauli 

Department of Theoretical Physics, Nuclear Physics Institute, Rez near 

Prague, CZ-25068 

E-mail: 



sauliQujf ■ cas . cz 



Abstract: Using the Green's function integral representation the Dyson-Schwinger equa- 
tions are solved directly in Minkowski space. Essential ideas of spectral techniques are 
discussed and applied on two renormalizable models: the Yukawa theory with massive 
pseudoscalar meson and conventional spinor QED. Within the momentum subtraction pro- 
cedure, the renormalization is performed analytically which leads to the usual dispersion 
formulation. The electron propagator obtained in this frame is compared with the solution 
of Euclidean Dyson-Schwinger equation and with the perturbation theory results as well. 
The proposed method has the advantage of obtaining solutions in both the space- and 
time- like regimes of momenta. In addition, when the coupling constant increases we find 
some unexpected discrepancy between the Euclidean and integral representation solutions. 
Especially for the supercritical couplings, the propagator pole is absent and there is no 
solution for spectral fermion function. 
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1. Introduction 

There are considerable interest in studies of strong coupling quantum field theories since 
many interesting phenomena, are believed have nonperturbative origins. In renormalizable 
field theory it is desirable to construct renormalization techniques which completely respect 
the symmetries of the underlying theory. In comparison with the conventional perturbation 
theory the nonperturbative treatment is more complicated and much effort is needed to 
perform renormalization procedure in a well-controlled manner. The fully regularization- 
independent method is proposed through the Dyson-Schwinger formalism, which 

is the main motivation of this study. The other motivation of the presented work is to 
avoid a constrains of Euclidean metric. Therefore the main part of our calculation is 
performed directly in Minkowski space and we offer the obtained numerical solution of the 
Dyson-Schwinger equations (DSEs) at all regimes of momenta. If the Euclidean analogue 
of solution is known and already present in the literature we perform the appropriate 
calculation of the Euclidean DSEs and compare with Minkowski solution. 
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Contemporary application of DSEs in hadronic physics offers rigorous insight into 
the infrared domain of Quantum Chromodynamics at zero temperature as well 

as DSEs provide solid foundation for nonperturbative continuous approach at nonzero 
temperature and density Q. Without the solution of DSEs it is difficult to reach any 
reasonable result on dynamical mass generation in the (Extended, Walking,..) Technicolor 
models 01,||8|,0]. In these models (usually without higgses) the values of the various 
condensates and the masses of the Standard model particle content are obtained from 
the solution of the appropriate gap equations. In practice, the system of the DSEs is 
truncated by an approximation of the Green functions that were thrown away. In addition, 
the truncated set of DSEs must be renormalized and as it is usually required by the 
method, the loop integral must be regularized before proper renormalization step. When 
any S-matrix element is completed from the Green functions the definite results must be 
l.renormgroup invariant, 2. gauge invariant 3. gauge fixing independent. All the three points 
are well understood in the perturbation theory treatment and it is desirable to build the 
similar nonperturbative method. The Pinch Technique |jT^,|ll],|12] and the Background- 
Field Method [ p!^ ] should offer the gauge fixing independent Green's functions, however 
how to fulfill all three aforementioned points is not answered by fully satisfactory way in 
DSEs treatment. In this paper, we particularly concern on the points 1 and 2 mentioned 
above. 

It is clear that the truncation of the DSEs system, if it is improperly performed, can 
violate gauge identity of the underlying gauge theory. In an Abelian versions of gauge 
field theory (scalar and spinor quantum electrodynamic) this problem was already solved 
by at least two technically different ways. Ball and Chiu [14| have derived the formula 
for the QED (truncated, proper) vertex that allows to close the DSEs system by a unique 
gauge covariant way (the lowest point Green's functions satisfy Ward-Takahashi identi- 
ties). The second successful approach is the longstanding method known as the 'Gauge 
Technique' |]l5[,[16|, |17], In this treatment the gauge covariant ansatze for the full 

(untruncated, i.e., not proper) vertex function is written in a terms of matter field spectral 
function. The technical advantage of the Gauge Technique is that the resultant equation 
for fermion propagator is linearized in spectral function. Due to this simplification, the 
Gauge Technique offers a solutions in compact analytical form. Of course, this linearization 
does not take place in the equation for the photon propagator and in fact this is one of 
the approach essential weakness: it cannot be true in general even for the case of electron 
propagator. Anticipate here, that no linearization does appear in this presented work and 
the appropriate integral equations are solved in their full form (they are third and second 
order integral equations). In the other side, the aforementioned linearization of Gauge 
Technique obtained equations do not exclude their reliability in the soft coupling regime 
and we shall mention the Gauge Technique once again when we discus infrared limit of 
electron propagator. 

In the last decade number of papers dealing with some improved vertices in QED 
DSEs studied the connection to dynamical mass generation |^],[^,|21], or they have 
been subjected to the various renormalization scheme used in its nonperturbative context 
|24|,p5[| (for a review of the earlier attempts see Q). In the paper |25], the modified 
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cut-off regularization metliod has been compared witfi tlie nonperturbative dimensional 
regularization scheme. This study exemplified the nonperturbative equivalence of different 
regularization-renormalization schemes. It also faces how careful the appropriate numerical 
procedure must be in order to obtain the reliable physical results (without breaking gauge 
and Poincare invariance of the underlying theory). We extend this studies to the case 
of (direct) momentum subtraction (MOM) scheme where neither numerical regularization 
is used. Likewise in the perturbation treatment, the all appropriate loop integrals are 
subtracted at certain value of external momenta and being then finite they are integrated 
analytically. Such a procedure leads to the usual dispersion relations for renormalized 
proper functions. We use the Green's function spectral representation which allows to 
convert momentum DSEs to the real equations for the spectral functions. At this point 



this approach is similar to the treatment used in the works |]22],|27| and of course, to 
the one used in Gauge Technique studies. Note at this place, that very similar method 
already has succeeded in the case of simple scalar models |2£],|3C] where the on-shell 
renormalization scheme was pronounced. 

It is worthwhile meaning that making some approximation the regularization indepen- 
dent and analytical answer can be obtained |25],|31],|^^. In very early stage of DSEs study 



there was made a sophisticated study [3^ of spinor QED where the propagators entering 
calculation are taken as bare ones and the resulting vertex turns out to be a hyper geometric 
function in that case. Furthermore, in rather special case (for massless fermion in rainbow 
and quenched Feynman-Fermi gauged QED) the authors of ||3^ found results analytically. 

In our paper the method of solution will be illustrated at two model cases. The first 
of them is the Yukawa Theory (YT), i.e., the theory of fermions interacting with spinless 
pseudoscalar boson. The Yukawa Lagrangian reads explicitly: 

- 90*^75^^ - /io^V> > (1-1) 

where mo , niQ^ are bare masses of fermion and meson corresponding to the unrenormalized 
fields f and $ respectively and go , ho represent the unrenormalized values of coupling 
constants. For the sake of simplicity, the Yukawa vertex is modeled by its tree value. In 
addition, it is assumed that h « g, therefore we also neglect the quartic mesons self- 
interaction. Remind here, that the quartic term in <I> is necessary due to the general 
requirement of renormalizability. The counter-term part has to cancel infinities appearing 
in the fermion loops contribution to scattering process — > $$. Here, the contribution 
from such process does not enter our DSEs due to their truncation mentioned above and 
an omission of ho^^ represents self-consistent approximation. 

The second model we employ is a conventional spinor QED. The approach is discussed 
for a general case of gauge covariant vertex. The numerical results are presented for the 
quenched, rainbow approximation leaving the complete treatment for the later presentation. 
Approximation used in this article has the advantage of being simple enough, which makes 
it an excellent testing ground for the proposed Minkowski analysis. The so-called ladder 
approximation (the bare vertex is used, the name follows from the ladder approximation 
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of the Bethe-Salpeter equation) is generally believed to be reliable in the Landau gauge 
only ^ = 0, and therefore we have used that gauge. Having all the numerical solutions 
stable and making comparison between Euclidean and Minkowski results we found that 
they agree only when the coupling is small enough. Increasing the coupling constant the 
obvious discrepancy does appear. Approaching QED coupling a = e^/(47r) to its critical 
value a = vr/S we are leaving only with Euclidean solution, whilst the spectral Minkowski 
equation tends to flaw. 

The necessary analytical formulas are reviewed in the next Section. It involves disper- 
sion relation technique, its connection with MOM and derivation of unitary equations for 
DSEs. The Section III is particularly devoted to the QED, where also the Euclidean version 
of fermion DSE is reviewed. Numerical results are presented in the following Section IV. 

Some details of the calculation and of the numerical method are explained in the 
Appendices. 



2. Spectral representation, Analyticity and Renormalization 

In the following sections we give overview of some basic facts about the Green's function 
Spectral Representation (SR), Dispersion Relation (DR) technique and their relation to 
the renormalization. 

The Lehmann representation [^],[^] for propagator can be derived from Lorentz co- 
variance and quantum mechanical requirement on the positivity of spectral density. The 
relations ( |2.1D display the necessary SRs for appropriate propagators entering the calcula- 
tion here: 



G{p') = / dcu ^ ' ' . = + dio- 

p'^ — LO + le — + le J p'^ — uj + le 

' — uj + ie p — m J — (J + 



p2 J p2 

GH/) = / dc^M^L = ^ + / (2,1) 

J p'^ — LO + le p'^ + le J p'^ — uj + le 

Here G, Sf are the full propagators of particles with spin and ^ respectively. G^^'^ describes 
propagation of massless particle which corresponds to the gauge field. The longitudinal ^- 
dependent part follows from the covariant linear gauge fixing qa = S{(,—df^A^) of quantum 
action. The single particle contributions ri6{uj — mf) are integrated out from the full 
spectrum a and it is assumed that remaining weight functions a's in ( |2.lD are smoothed 
ones and not complicated distributions at all. 

It is notable in this place, that there is less formal proof of SRs for the simple Quantum 
field models. These models are represented by the scalar quantum field theories without 
derivative interaction. The existence of integral representation was proved to the all orders 
of perturbation theory even for an arbitrary n— points Green's function. In this case, the 
so called Perturbation Theory Integral Representation (PTIR) was derived by Nakanishi 
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|37]. Furthermore, it was shown that the PTIR is unique, which property appears to be 
very useful in some apphcations (see for instance |28],|3C], and references herein). 

Although the formal derivation of Lehmann representation is rather straightforward 



(see some standard textbooks |35], ||3^), it necessarily breaks down when considering a 
theory without free particle asymptotic states, i.e. the theory with confinement. If the 
confinement takes place in a given theory then the particle can never be on mass shell 
and the appropriate propagator should not posses singularities at the real time-like axis 
of p^. The absence of Lehmann representation should be a good signal for confinement 
pO| ] , 1^ . Actually, as we will show in the Section devoted to the strong coupling quenched 
QED, the disappearance of solution for Lehmann function rather sharply coincides with 
the transition to a confined phase. 

As was mentioned, the nonperturbative regularization - renormalization procedure is 
not so transparent as in the case of perturbative treatment. Our task is whether the 
renormalization procedure should be performed analytically in an easy and transparent 
fashion as it is in the perturbative approach. First of all we must answer the question that 
naturally arises: what are the physical criteria, which will one to determined the solution 
of DSEs as a physically meaningful. In particular, the following is required for two point 
function: 



(1) The renormalized full Green's function satisfies its own renormgroup equation 7 = 
/i^G(p^,/i), where gamma represents logarithmic differentiation of conventionally defined 
field strength renormalization constant Z (explicitly introduced later by Eqs. (p.3|)) 7 = 
— yU^, i.e., the unrenormalized propagators 5/o, Go, Gq^ have to be manifestly independent 
on the choice of renormalization scale. 



(2) The on-shell renormalization scheme (ORS) should be a special choice of the renor- 
malization scale fi = m. Furthermore, the position of the pole p'^ = m? is renormgroup 
invariant quantity. 



(3) The Green's functions are an analytical ones. Up to a subtracting polynom, the 
real and imaginary parts of proper Green's function (one particle irreducible diagram with 
truncated legs) is uniquely defined by a dispersion relation. 



(4) The off-shell proper Green's functions calculated within DSE approach should 
admit renormalization. It implies that the possible maximum number of subtraction is two 
for self-energy and at most one for the possible triplet and quartic vertex renormalization. 
All the coefficients in the appropriate subtracting polynom should be absorb-able in the 
counter-term part of the original Lagrangian. 



(5) As far as it is possible, the renormalization procedure should respect classical sym- 
metry of the theory. Particularly, when dealing with gauge theory, the Green's functions 
should satisfy Ward-Takahashi identities. 
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In the next subsections we review the general framework of analytical renormalization 
technique and explain the renormalization scheme which is actually demonstrating that 
the Green's function obtained by this technique satisfy all the above requirement 

2.1 Renormalization of DSE for Yukawa pseudoscalar 

Here we start with the discussion of boson propagator. The unrenormalized version of DSE 
reads 



GoHp^) =P^-ml-Uo{p^) 
no(p') = iglj Tr^TYoip,l)So{lh'Soil-p), (2.2) 

where the index represents bare quantities. To renormalize YT we conventionally in- 
troduce the renormalization functions with the appropriate counter-terms. Choosing some 
arbitrary renormalization scale fi'^ they are: 

m^o = Zm^m^ifi); 6^^ = ml - m{fif , uiq = Zmm{^i); 6m = mo - m{n) 

^1/2 ^1/2 

90 = -^9 ,ho = ^h , (2.3) 

where Z$ (Z,j,)is a renormalization boson (fermion) field-strength constant , Z^, is the 
mass renormalization constant, m^ is a bare mass while m(/i) represents renormalized 
mass. Furthermore, Zh and Zg represent the renormalization constant of quartic <I>^ and 
triplet vertex respectively, i.e. for instance we can write for the renormalized proper 

triplet vertexTy (p,p — k) = ZgYyoiPiP ~ k). First of all we renormalize the propagators 
with respect to the field renormalization only i.e. Go = Z(f,G, Sq = Z^S. For this purpose 
we multiply the unrenormalized DSE (^]^) by the constant Z^. A simple algebra gives 

G'Hp^) = Z^(p2-mg)-no(/) 

mp') = ig'l Tr^Ty{p,l)Sf{l)^^Sf{l-p) , (2.4) 

where the proper vertex Fy satisfies its own DSE. To truncate the system of DSE we 
make most simple approximation Fy(p, /) = 7^. Since the pseudoscalar particle $ requires 
quadratically divergent mass renormalization and logarithmically divergent field strength 
renormalization the relation between renormalized and unrenormalized self-energy function 
must be of the form 

no(/) = no(M') + n[,(M')(/ _ + n(^;/) , (2.5) 

where the renormalized self-energy satisfies double subtracted DR 



n(/u;p') = j du 



(jp' — oj + ie) 



P - 
a; — /i^ 



(2.6) 
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Of course, from the Rel. (|2.5| ) we obtain the standard receipt for calculation of the propa- 
gator in MOM scheme 



dp 



n(//;/) = no(/) -no(//2) - 42^oiP%2^^2(j>^ - fi^) . (2.7) 



From Eq.(^.5|), one can readily see that the unrenormalized value no(^^) corresponds to the 
dominant part of the mass renormalization constant and the derivative nQ(//^) corresponds 
to the field-strength renormalization, explicitly we have 

Z^ = l + 



Note that in the perturbative context this schemes sometimes referred as the BPHZ 
renormalization scheme and the appropriate subtraction procedure is called (Bogoliubov) 
R-operation. In the next text we drop out the labeling MOM since no other renormalization 
scheme is used throughout this article. 

Employing the well known functional identity for distributions 

^ P.—^ iTx6{x-x) (2.9) 



we see that the function 7rp{p'^) represents the absorptive (imaginary) part of renormalized 
self-energy n(p^) as well as of the unrenormalized one. The last statement is unrelated with 
the question of (in)finiteness of the counter-terms, since such renormalization procedure 
has nothing directly with the presence of infinities and can be consistently applied to a 
theories which are UV convergent. The detailed derivation of the weight function is 
relegated to the Appendix A. Here we simply review the result: 



(JLf 



00 

1 



+mrf J dp [-{oj -m'' - p)a^{P) +mas{l3)] X{uj;m^,P) 



00 00 

+ I da I dp (^^ 



^{oj - a - P)a^{a)a^{p) + as{a)as{P) ) X(u;; a, /3) 1(2.10) 



(m+m^)^ (m+m^)2 



where the function X can be expressed through the triangle function A by the following 
way 

X(a; b, c) = ^^^!(^e(a - (^/(6) + v^c))^) (2.11) 
and m is a pole mass of fermion. 
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A boson pole mass rrifj, = m^{m(f)) is conventionally defined as G "^{rn^) = 0. From its 
definition it reads: 

. ■ r 2 2 



00 — m?. 



(2.12) 



Due to the algebraic simplicity of the DRs in MOM scheme we can immediately rec- 
ognize that the inverse propagators G renormalized at two different scales ^' just differ 
by finite polynom a fin. + b/mP^- Written explicitly, the relation reads 

" ~ J {U - /i'2)(a; - ^2)2 

bfin. = m [1^) - m [n) ^ 5 

= m{fi)-m{iJ,)+ du;-^ ^- (2.13) 

The MOM identities (??) exhibit the evolution of the propagator within the change of the 
renormalization scale. 

2.2 Renormalization of DSE for Yukawa fermion 

The extension to the fermion case proceeds similar way. DSE for fermion propagator reads 

57i(|^) = Z^[^-mo]-So(p), (2.14) 

where ttiq represents bare mass of fermion and the renormalization constant was already 
introduced in ( [2.31 ). It is convenient to split unrenormalized self-energy to its dirac vector 
and dirac scalar part 

So(p) =^ao(p2)+6o(p2). (2.15) 

The scalar functions oq, ho can be easily identified from the explicit expression for fermion 
self-energy 

So(i^) = -ig' I (^ry (p, i)Sj{ih'G{{p - If) . (2.16) 

The fermion self-energy is only logarithmically divergent and one subtraction is sufficient 
to make the scalar functions oq, bo finite. The renormalized self-energy then reads 

^ifi,p)=Mf^,P^) + Kf^,P^). (2.17) 
where the subtraction leads to the following DR's for renormalized a, b 

2^ _ „ /^2^ ./„2^_/",„ Pv{s){p^ - fl^) 



a{fj,,p ) = ao{p ) - ao(/U ) = J ds 
6(/x,/) = 6o(p')-6o(/^')= / ds 



(p2 — s + ie){s — /x2) 

p,(s)(p2 -^2) 



(p2 — s + ie)(s — /i2^ 



(2.18) 
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Then the renormahzed version of ( 2.14| ) can be written hke 



SJ\p) = A{^^)[J^-m{^I)]-J:{^I,p) . (2.19) 

where m(/x), represents renormahzed fermion mass (coefficient of /^) fixed at the 

scale 11. As it is usuaUy in DSE treatment the relation ( 2.19| ) can be equivalently rewritten 
into the form: 



A(p'^) = A(;u) -a(/i,p2), 



(2.20) 



where we do not indicate explicit dependence on /i in the renormahzed functions S, A, B 
for brevity. For the absorptive parts 9^ = 9a = -Kp^^'^B = 96 = ups we can find the 
following results: 



Pi,(w) 
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(47r)2 



+ '^f J dPp^{P)Xi{u;;P,m ) + J dP J dap^{(3)a^{a)Xi{uj; P,a) 

4m^ 4m^ (m+m^)2 



(2.21) 



(4^)2 



00 

r^rfmX{uj;ml,m'^) J daas{a)X{uj; 171^,0) 



+ 



J d(5p^{[i)X{uj-f3,m')+ J dp J dap^{P)a,{a)X{uj-a,P) 



, (2.22) 



where the function X has been introduced earlier and Xi{x;y,z) = X^^'^{x;y, z){x — y + 
z)/{2x^)e{x - (v^y) + x/(z))^). The detailed derivation of Eq.'s (l2l2ll), (pT2^) is presented 
in the Appendix A. 

The renormgroup invariant mass function which is conventionally defined as 



Mip^) = B{p^)/Aip'). 



(2.23) 



Since the on shell mass m = m{m) is usually the one which is the best known experimen- 
tally, it is useful to introduce the relation between m and the one renormahzed at p. From 
its definition S{m)^^ = we have 



We explicitly choose 



m = A{p) [m{p) + a{p] m )m + b{p; m )]. 



A{p) = 1 



(2.24) 
(2.25) 
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in this paper and Yukawa fermion mass m(^) is fixed through ( 2.24| ) such that m = 1. 

Imposing the renormahzation condition on S at two different scales fi, and /i' we can 
again recognize that the inverse of the fermion propagator differs by certain finite piece 
c/m.+ jjpdfin.-, here we only note that the coefficients Cfin.,dfin. can be expressed through 
the functions a, b{fi; /i'^) and a, b{fj,'; fj?). 

2.3 DSE for a's- The unitary equations 

The derivation of the DSE for Lehmann weights is presented in this section. Evaluating 
the imaginary part of the trivial identity 1 = G{p'^)G{p'^)~^ with SR (2J.) used for G and 
DR used for 11 in the DSE, i.e. in eq. G~^ = p^ — m^{p) — Il{p;p'^) yields 



a{uj){uj -ml{p)) = 



p4^) 



ijj — nij 



+ {p<f>*cr}, 



(2.26) 



where the symbol {p^ * a} represents following principal value integral: 

p{s)a{x) + a{s)p{x)^^^ 



P. / dx- 



The pseudoscalar propagator residuum 



lim 



s — X 



1 1 

p - "li 



G ^(p2) 

is most easily evaluated through the appropriate DR for self-energy 



(2.27) 



(2.28) 



1 - n'fmr 



1 



uj)(m\- p^){m\^ p^ 
(m| — ijoYiuj — p'^Y 



(2.29) 



It is not surprising that the equation ( |2.26D looks particularly simply in on mass-shell 
renormahzation scheme: 



(JiUJ 



(a; — oo — m\ 



2 ■ 



(2.30) 



The fermion case can be treat by a very similar way. From the fermion DSE we can obtain 
for the residuum 



lim ^_ = [1 — o'(/i; rn?)m — b'{p; rn?) — a{p; m^)] ^. 



(2.31) 



where we explicitly used the renormahzation condition ( 2.25 ). Furthermore, writing the 
trivial identity S~^S — 1 = in a suitable form: 



m 



+ / duj- 



m{p) 



p"^ — uj + ie 

{^Pv{s) + Ps{s)){p^ - p^) 



ds 



{p^ — s + ie){s — p^) 



(2.32) 
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then we can easily arrived at the relations between Lehmann weights cr's and the absorptive 
parts of self-energy p. Projecting the obtained result by and ^J^^i t> leads to the 
coupled set of integral equations 

I ^ fi+m{p)f2 , , m{p)fi+ujf2 
<7t,(a;j = j^—- ; crs[uj) = 

_ ujpyjuj) +mps{uj) 

h = r 5 >rUj[(7v* Pv\ + [(7s * Ps\ 

_ mpy{uj) + ps{uj) 

h = T o V * Pv\ + \<7v * Ps\ , (2.33) 

LO — 

where w is (positive in our metric) time-like momentum u = and where we have used 
the abbreviation for the real functional: 

2 



oo ^ p{s)cj{x) +a{s)p{x)^ 
s — X 



[a*p] = P. I dx'—— — _ (2.34) 



where P. stands for principal value integration. 

Note also, that the fixed renormalization implies the unique determination of the prop- 
agator residuum. Henceforth, its value ca be obtain obtainable without on-shell differenti- 
ation of the self-energy function. Clearly putting p t any fixed value the residuum can be 
easily extracted ( simplifying choice is for instance ^ = or = 0). 

3. DSEs in Quantum Electrodynamic 

After the more general introduction we describe Minkowski formalism which is necessary 
in QED DSEs treatment. For this purpose we used the conventions already established in 
the last two previous sections, the differences that appear are emphasized. In the end of 
this section we also review the Euclidean fermion DSE, its solution serves us for numerical 
comparison. 

In general, there is no isolated pole but propagator singularity coincides with the 
branching point. The analytical structure of QED fermion propagator was the subject of 



the initial study ||4^. The authors of the paper [^] converted the integral gap equation (for 



its explicit form see Rel. ( 3.15 ) in the text bellow) to the non- linear differential equation 
which has been then solved by graphical method. In addition, they proceed the backward 
Wick-rotation of their equation to the time-like regime of momenta. For a large coupling 
enough, they did not find zero in the inverse electron propagator (for a real p^). This 
unexpected disappearance of physical branch point leads the author to the conclusion that 
confinement should exists even in the Abelian gauge theory. Having the ultraviolet cut- 
off A fixed to the some finite but large value (A >> M(0)), the authors of |4^] identify 
the coupling of the phase transition to be exactly the one related to the dynamical mass 
generation, i.e. a = 7r/3 in Landau gauge. 

In fact, the claim that the Abelian gauge theory has a confinement phase sounds 
strange from the conventional wisdom based on the continuum Abelian gauge theory. It is 
clear enough, that the electrons must be a free particles when the coupling is small enough. 
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The most recent work can divert this unrehabihty. It was argued in the paper [^] that 
the Abehan confinement should exist only at the above certain value of strong coupling 
constant, probably where the photon reveals the dynamical mass. Making the quenched 
approximation, the photon mass generation can not be answered, but the corresponding 
disappearance of the physical branch point is observed and confirmed in our study. In 
agreement with our expectation, this transition is observed for the value a ~ 1 of QED 
coupling constant which coincides (up to the numerical accuracy) with the failure of our 
spectral approach. 

The later analyze ||4^ founded that the electron propagator (in the ladder approxi- 
mation) develops singularity for the all couplings. The results of the paper [45| support 



partially the conclusion made in ||4J] (but, not for all the couplings): The electron propa- 
gator has not only one real branch point, as it is physically expected, but it also embodies 
two other complex conjugated singularities (the position of them entails troubles with the 
analytical continuation, see the discussion in (P). To conclude, note that it is generally 
believed, that this unexpected pole's complexity should vanishes for the exact solution 
(particularly when the coupling is small enough). In fact, the existence of the complex 
branch points is questionable and the appropriate answer depends on the approximation 
employed. Note that, it is not necessary in the contradiction with the result of us and with 



the earlier study of Fukuda-Kugo [42|. It would be only in the case if the pure real pole 
solution does not vanish for the coupling constants large enough. 

In the spectral approach used thorough this work, the dominant part of the propagator 
is driven by the real pole part: rj(p^ — m?)~^ and the changes coming from the interaction 
are involved in continuous part of Lehmann spectrum 

~ u! + ie 

which may but need not to be analytical at the branch point p'^ = m^. Recall at this place 
the Gauge Technique (zeroth order iteration) result [^] where in class of covariant gauges 
the infrared behavior reads 

• (3-2) 

We actually see that the interaction partially suppressed pole singularity for those gauge 
fixing parameters that are less then Yennie gauge (^ = 3 in Yennie gauge) . 

Furthermore, in contrast to the YT , we are now dealing with the gauge theory and the 
appropriate renormalization scheme should respect the gauge identity. Mainly the vacuum 
polarization tensor should be transverse 

'?^n^.(g) = (3.3) 

and the proper photon-electron-electron vertex should satisfy Ward-Takahashi identity 

S-\p)-S'Hl) = ip-l)^T^^{p,l), (3.4) 
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which uniquely determined the longitudinal part of the vertex |14] 



+ 



2 p'-P 
Using the DRs for the function A,B 



p^ 



'B{p') - B{f)) . (3.5) 



1 



da 



1+ da 



{p^ — a + ie){a — /i'^) 



(p"^ — a + ie){a — /i^) 



(3.6) 



we obtain the integral representation for 



Ti{p,l)=Y I 



duj- 



u; — /i^ 



{p^+in{p+ /) 



(pM + ZM) j 



duj 



dio^ 

p"^ — UJ 



duo 



P-UJ 



dio- 



(p2 



to 



(3.7) 



(p2-a;)(Z2 

From this expression one can read immediately see that the longitudinal part of QED 
vertex is free of any kinematical singularities and we also see that the only coefficient 
of is explicitly renormalization point dependent. Furthermore, we can note here that 
this should be the whole fi explicit dependence of the full vertex F^ = F^ + F^, since its 
transverse part (satisfying {p — 1).Tt = 0) must be finite (for more details of (perturbative) 
MOM renormalization scheme used in QED see |]46|| ). 

Substituting the gauge covariant vertex (^^) into the DSEs for fermion and photon 
propagator we obtain two closed equations that can be solved numerically after the renor- 
malization (see most recent paper |^l| and references therein) Adopting rather standard 
notation: let the constants Zi ,Z2 and represent the vertex, fermion wave function and 
photon wave function respectively then the renormalization proceeds by the standard way: 
the Dirac and scalar function ao(//) and 6o(m) iii unrenormalized Eg must be absorbed in 
Z2 and Zm renormalization constant by the same way as it is happen in the case of Yukawa 
fermion propagator. Furthermore the unrenormalized vacuum polarization n(/u) should be 
absorbed in the renormalization constant Z^ and the infinity of F'^ is canceled against the 
constant \ — Z\. The renormalization scale used to renormalize vertex and the electron 
propagator should be the same due to the WTI ( |3.4D that requires Z\ = Z2 { Again we do 
not explicitly indicate the dependence on the renormalization scales (i.e., Zi = Zi{p))). 

In our MOM 'imaginary part analysis' we must self-consistently obtain DR ( p.6| ) for 
fermion function. Further we should obtain DR for renormalized photon vacuum polariza- 
tion 
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oo 

r 2 2 

MOM ( ,,2. J2\ _ / J, , 9 — 







where its absorptive part reads 



7rp((j) = + 2mV^^)Vl-4m2/a;e(tJ - W) + O(e^) (3.9) 

with aQED = e^/(47r), e is a charge of electron, m is on sheh electron mass. Without 
making some other approximation the full treatment with the DSEs leads to the evaluation 
of large number loop integrals and it also requires a careful numerical treatment due to 
the presence of Landau ghost. The complete solution of this problem is relegated to the 
forthcoming paper [47|. 

Using the bare vertex and quenched approximation the fermion DSE reads 

Z,[^-mo]-ie' I ^G^,''ik)rS{p-k)Y = 0. (3.10) 

and straightforward calculation gives a{fi;p'^) = 0. This entails equality A(fi,p'^) = 1 for 
all square of momenta |48] if the condition A{fi) = 1 is imposed. The unitary equations 
are notable simplified in this case: 

. X fi+m{fi)f2 . m{fi)fi+ujf2 

crv[uJ) = ' ^H^J = Of \ ' 

UJ — m'^(p) w — m^(/uj 

_ mps{uj) ps{uj) 

fi = r 2 + * Ps] ; /2 = r 2 + Ft; * Ps] , (3-11) 

to — UJ — 

where the absorptive part of the renormalized self-energy 9S(p) = h{p;p'^) = irpsip'^) is 
given by 

- • (3.12) 

Some details of derivation of Rel. (|3.12| ) are given in the part b) of the Appendix A. 

In order to make a careful and constructive comparison between Minkowski approach 
presented in this work and the standard Euclidean formulation we should compare with 
some known Euclidean results presented in the literature [^],[ 22 1. In space-like region the 
Eq. ( p.lO|) is transfered to 



pM = -3 (l - ^) + £^ daa,i<^) (l 



B{x) = mo + ^ j dyK{x,y) 



B{y) 







Kix,y) = \ (3.13) 

x + y + ^/{x-yy 

where Wick rotation and angle integration have been done and where we have used Z2 = 1, 
X = p\ = —p'^,y = k'^ = —k"^- The renormalized equation then reads 

00 

3« l\ ^r<, ^ ^(C>y) 



B{C,,x) = m(C) + ^ y dyV{C,,x,y) 



V{C,x,y) = K{x,y)-K{C,y) (3.14) 
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where C is square of renormalization scale —n^- Choosing the scale C to be zero and further 
scaling the renormalized m{C) mass as m(0) = 1 leads to the particularly simple expression 
for (|3ll| ) 



m X) = 1 + ^ / - l) -^g^ (3.15) 







(for 'massive photon' case see Eq. (3.14) in the work |42|). Stressed here, that due to the 



Landau gauge the solution of Eq.( p.l4 ) is only slightly deviating from the solution with 



Curtis- Pennington vertex implemented ||2^. From the paper |25] the asymptotic behavior 
of dynamical mass is known analytically. Henceforth, any reliable solution of electron DSE 
must behaves at ultraviolet like 

M(p|) = M(C)(^) . (3.16) 

4. Numerical Solution and Results 
4.1 Yukawa theory 

The resulting coupled nonlinear integral equations ( 2.33| ),( p.26 ) for the functions a^^s and 



£70 require the knowledge about the value of physical masses ( 2.24 ), ( 2.12 ), the propagators 



residua ( |2.29| ),(??) and the complete knowledge of the absorptive parts of self-energies 



( 2.10 ),( 2.21 ),( py2^ ). The equations have been solved by the method of numerical iteration 
which seems to be particularly useful for this purpose. The one loop perturbation theory 
result is used as the zeroth order of this iteration. Then, several hundred of iteration steps 
have been proceed to achieve a stable solution. For purpose of numerical integration we 
choose the Gaussian quadrature method. Taking a reasonable number of the integration 
mash points and adopting the principal value integration described in the Appendix B, 
then the whole numerical procedure is rather stable against the changing of mesh points 
density and their number as well as. 

The all presented results for YT are evaluated for the zero value of renormalization 
scale /x^ = which choice is common for propagators corresponding with both the particle 
of YT. The physical mass is usually the best known experimentally and we face that this is 
our case. Let us assume that experimentalists (living in our toy model world) found their 
values: = 0.15m, m = 1. In this case, we are enforced to calculate m(ij,) and m^{^) 



from the Rels. ( |2.24| ), ( p. 12 ) which procedure does not cause any troubles in our iteration 



method. The above described procedure has clear numerical advantage: the branching 
points lie at the values of momenta which are fixed at each iteration step, but it has also 
disadvantage: the comparison with some results obtained in the Euclidean formalism is 
not so straightforward (note, m(;u) is fixed for some space-like /i^). To make the most 
accurate comparison with Euclidean result, we should fix the renormalization scale and the 
renormalized masses at the same values that are used in the Euclidean DSEs. We prefer 
the first scheme in YT since the appropriate Euclidean solution is not published elsewhere 
(but we use the second approach in QED case). 
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Fig.l shows the absorptive parts of self-energies, i.e., the functions Pv,Ps for fermion 
and Pfj) for pseudoscalar. They are plotted for the following values: A = 0.1,0.2,0.3 of the 

2 

coupling strength defined as A = For a large momenta the spectral function P(f,{x) 
grows linearly with the square of momenta which is a consequence of quadratic diver- 
gence in Ho. Due to this, what we actually plot is the rescaled function p^(x)/x. All 
the functions start to be non zero from their perturbative thresholds- QA{s),QB(s) from 
s = {m(j, + m)^ and 9n(s) from s = Am?. Because the negative parity of the field ^> 
the functions Ps-, Ptf, (and (7s,cr^) are positive while p^ (a"„) is negative. The quadratic 
divergence of the unrenormalized self-energy 11 leads to the quadratic dependence (-|-cor- 
rections) of the renormalized 11. (in fact the recent models (see for instance [^],[^]) of 
particle interactions attempted to avoid of quadratic divergences that necessarily follows 
from Yukawa sector of Standard model) The quantity Q^/^ of dimension [mass]^ defined 
with the help of eq. Q{p^) = 'm'^{p) + Il{p;p'^) can not be called the dynamical boson 
mass, since Q becomes negative above the certain value of Euclidean momenta. The same 
happens at the time-like axis for ^Q{p'^) (while QQ{p'^) = Trp^ is positive as it is clear from 
Fig.l). The appropriate behavior of the function Q is displayed in Fig. 2. The momentum 
dependence of the fermion functions A, B is dominated by the (perturbative) logarithm of 
p^. Indeed, the Yukawa fermion propagator is almost given by its free form corrected by 
small perturbation. Fig. 3 and Fig. 4 display the momentum dependence of the function 
A, B for space-like momenta and time-like momenta respectively. Although, the smallness 
of the ratio mphi/m is motivated by the desired enhancement of fermion self-energy and 
subsequent suppression of pseudoscalar one, nevertheless, we still see that the fermion- 
meson loop becomes irrelevant perturbative contribution for all studied couplings. We 
display the dynamical mass M = B/A in the infrared domain in the Fig. 5. since this is 
the mainly interesting regime of momenta. The appearance of the critical coupling Ac is 
the consequence of the quadratic momentum dependence of pseudoscalar self-energy. Its 
value slightly depends on the numerical cut-off A used in our unitary equation, noting that 
within our numeric (A^ = lO'^m^) the solution of unitary equations fails at Ac — 0.3. 

4.2 QED fermion propagator 

In contrast to the model discussed previously, the strong coupling QED is less driven by 
the perturbation theory (Recall here the famous paper [^] dealing with scalar electrody- 
namics). Before presentation of our numerical solution, let us recollect the main results: 
1. Comparing the electron propagator obtained from the solution of the Unitary Equa- 



tions (UEs) (3.11) with the propagator calculated in the Euclidean formalism we find that 
they never exactly agree. The exception is the case of very small coupling a ~ 10~^, where 
both approaches seem to be equivalent (up to the numerical errors). Stressed here, that in 
this case, they are almost indistinguishable form the perturbation theory result. 

2. Previous statement is valid also in the case when the small photon mass A is intro- 
duced. The obtained results are then slightly changed quantitatively and the disagreement 
discussed in 1. is somehow small (the Euclidean x Minkowski results numerically agree, only 
when A ~ m, which is the case, we are not interested in). It leads us to the conclusion, 
that this discrepancy do not fully follows from the masslessness of the photon and from the 
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appropriate subsequent coincidence of the electron propagator singularity with the branch 
point. Note, that the last property is not simply fulfilled for nonzero A. 

3. The Minkowski (spectral) solution exists only for subcritical regime of aqED, whilst, 
as it is well known the Euclidean solution can be easily find for both the sub- and super- 
critical couplings. To be more precise, we have not found any solution of UEs when the 
coupling was larger then Ocs. = 0.95 (c.s.=critical spectral). Further, making an estimate 
of the momentum DSE solution (not assuming spectral decomposition) we have found 
strong evidence for confinement for the coupling larger then a = 0.91 it 0.1. (Remind, that 
the critical coupling in QED is defined such that B{p'^) = for mo = 0, (chiral symmetric 
phase) for a < Oc and the solution for B{p'^) is non-zero and finite for a> ac', mo = when 
momentum cut-off A is implemented (without finite A the function B tends to diverge). 
Its value is known from the quenched, rainbow study, where ac = , more sophisticated 
studies do not deviate significantly from this value ). 



The UEs (|3.1lD complemented by the equations for the absorptive part ps and for 
the residuum have been solved by the method of iterations. In order to have an infrared 
behavior under the control and in order to see the effect of the vector boson mass we also 
implement the small mass parameter A into the photon propagator: 

. Restricting to the small photon mass case, the Dirac part coefficient function A is 
approximated as A{p'^) = 1, i.e. as it would be in the massless case. The expression for 
function ps is then slightly change (see Appendix A), noting here that the limit A — > 



is smoothly achieved as it follows from the appropriate relations (B.7) presented in the 
Appendix A. The solutions with several different photon masses have been obtained. Then, 
the dynamical mass M = B at space-like momenta is calculated from the absorptive 
functions ps,crs,(T^, which are primary solutions here. We can see from the Fig. 6. that 
the effect of the small photon mass is almost irrelevant even if the coupling is relatively 
large (a = 0.6). The space-like renormalization scale is chosen to be large //^ = —10^ 
when compared with renormalized mass rn?'{p) = 400^^, which choice corresponds with the 
Euclidean solution already presented in the work . The appropriate Euclidean solution 
of Eq. ( 3.14] ) is added for the comparison. The same (in) dependence on the parameter A 



is observed for any choice of p,m(p). In order to achieve good numerical stability of UEs 
some very small photon mass is always used in. The presented results in this work are 
calculated with A = lO^^m. 

Solutions were obtained for the Fukuda-Kugo equation (FKE) ( |3.15| ) for the couplings 
from 0.01 to 2.0. The expected damping p2| , ^] of dynamical mass to its negative values 
was observed for supercritical couplings a > 7r/3. Using the same renormalization choice 
M(0) = 1 the UEs were solved. The resulting solution for M is plotted for the following 
couplings: 0.2,0.4,0.6 and compared with the FKE solution. All these solutions are plot- 
ted against the space-like momentum in Fig. 7. From this we can see the relatively large 
discrepancy between the results of FKE and UEs. Making great numerical effort we have 
found the main cause numerically. Let us assume that the residuum and the principal value 
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integral integrations that appear in UE are somehow over-estimated due to the infrared 
enhancement of the Lehmann weights. Let us introduce small coupling constant dependent 
infrared cutoff c * to the UEs: 



h ^ r^^^ + K * p.Mco - m'il + c)) 



/2 = + [a, * psMu - m\l + c)) 



UJ — m 

2 



m{ix)m — ij? J \ J UJ — p? 



a 
2^ 



(4.2) 



Note that taking c = then the modified unitary equations (MUEs) (^^) clearly reproduce 
the originally UEs. Looking at the solutions of MUEs we see that their solutions are in the 
reasonable agreement with the solutions of FKE. Fig. 8 displays the mass behavior at the 
time-like regime obtained for the same choice of renormalization mass as previous. The 
maximum corresponds correspond with the mass-shell point . Note that their existence is 
in contradiction with the assumption of confinement |^^. Although, we are pretty sure 
that the solution must look like as in this figure, in order to assure, we solve the time-like 
continued equation of Fukuda-Kugo. We were not able to find the exact and fully complex 
solution numerically, however, we made a sophisticated estimate of a real part of B{p'^) by 
the principal value integration. It is achieved by the standard numerical introduction of a 
finite epsilon, here it is chosen to be a fraction of pole mass and presented solutions are 
obtained with the choice e = 0.03m^, i.e., the singular integrand is approximated as 

" ' . (4.3) 



'z-B^z) {z-B^{z)y + e 

where z is a time- like momentum z = p^. 

The soft coupling results are added to the Fig.8 and fully agree with our expectation: 
the branch point always exists because the mass function B'^{x) always cuts the graph of the 
function y{x) = x. This statement is valid for QED with subcritical couplings a < 0.915. 
At most, we expect about ten percentage deviation from the exact value of acrit which 
should be a consequence of the numerical method weakness. Furthermore, we can see 
from the Fig.8 that FKE solution for ^B reasonably agree with the solutions of MUEs, 
especially when the coupling is small enough. But we have rather different situation in 
strong coupling QED. The maximum of the mass function disappear and the mass function 
never cuts the graph y/x. The supercritical coupling quenched QED is a confining theory, 
the electron propagator has not a branch point and is free of any singularities. In this case, 
e introduced above can be safely limited to the zero value and we actually omitted it when 
FKE was solved for a > acrit- The results are plotted in Fig. 9 where we also add some 
subcritical solutions for better comparison. The subcritical solutions have been obtained 
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with the help of time-hke continued FKE and there is no spectral solutions presented in 
this figure (the value of infrared cut-off c makes the MUEs untrusted). 

The program using the Euclidean equation ( |3.14| ) was run with several coupling a from 
0.2 to 8.0 (!), and with the renormalization choice m(— 100) = 10. Whenever a coincides 
with the one published in |22], then the results here and the ones obtained in |22| should 
be numerically identical. Again, we see for supercritical couplings, there is a small region 
where the dynamical mass is negative. When the coupling constant is extremely large (say 
a = 8) the negative damping becomes a relatively fast oscillation around the zero mass 
axis. This feature should be in agreement with |22] but in a subtle disagreement with the 
paper [^], where the oscillation appears to be purely positive. The solution of B in space- 
like regime is presented in Fig. 10. The solutions of MUEs and of the Euclidean momentum 
DSE are compared only in the space-like regime. The MUEs solutions are presented in in 
Fig. 11 for time- like momenta. 

An unexpected feature of the spectral equations is that for any set {a > 0.9, /x^, m(//)} 
the absorptive self-energy becomes oscillating around the zero axis. This behavior was 
verified to be insensitive to number of grid points, but was not stable against the infrared 
cut-off introduced above. In fact, the positions of minima's and maxima's are slowly walk- 
ing when the cut-off c is varied by hand. The significance of this for spectral representation 
and its connection to QED is not completely understood. One possibility is that it may 
be signal the failure of Lehmann representation. The examples of this feature is shown in 
Fig. 12. It should be emphasized here, that all the results for M(p^) obtained by solving 
of MUEs and presented in this work have been calculated from non-oscillating, smoothed 
(soft coupling constant) p^'s. 



5. Conclusions and outlook 

In this article we convert the Dyson-Schwinger equations for propagators to the real equa- 
tions for their Lehmann weights. The possibility to do this, is based on the assumption of 
the existence of Green's function integral representation which should be valid in nonpertur- 
bative regime too. They renormalized equations have been solved without some unwanted 
linearization or angle approximation. This is the novelty and significant advantage in com- 
parison with the usual manner commonly used in the literature. The developed approach 
has been tested on two models. The QED electron propagator has been compared with 
the results obtained in conventional Euclidean formalism. We found certain discrepancy 
between these approaches which appears to be rather large when the coupling approaches 
its critical value. Furthermore, the solution obtained within the help of Lehmann represen- 
tation is fully absent for supercritical couplings regime of QED. In this case the physical 
pole propagator singularity disappear and we argue that ladder QED does not describe 
free electrons (positrons) at all. 
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A. Evaluation of the off-shell dispersion relations 

In this Appendix we derive the appropriate DR's for the proper Green's functions of YT. 
As a first step we analyze the fermion self-energy. 

A.l Yukawa fermion self-energy 



As follows from ( p. 16 ) the general expression which has to be considered has the following 
structure: 



(2^)4 (/2 _ a + ie){{l - p)2 _ /3 + ie) 

= -7^ j'dx f -C.^(l-x)+a 

^ J J i2Tr)^P+p^xil-x)-ax-Pil-x) + i€]'^ ' ^ 



where the whole prefactors are absorbed into the formal symbols Ca^b ■ Written in the 
terms of continuous functions a^^s and Lehmann weight a for pseudoscalar, the symbol Ca 
and Cb are identified as 

Ca = J da {6{a - m^) + a^{a)} J dp {6{P - m^) + a{p)} 

Cb = j da [a^'^5{a - m2) + a,(a)} j dp {5{P - m^) + a{p)] (A.2) 



Noting that they should be collocated in front of the fractions in (A.l). 

Observing the Dirac structure, the unrenormalized functions ag, ho can be simply iden- 
tified. After making subtractions ( p. 18 ) we can arrive at their renormalized forms 



a{fJ.]P^) = ig^Ca I dxdy 



dH (1 - xfx2{p^ - fi^) 







KfJ-iP^) = -ig^Cb / dxdy 



(27r)4 [/2 + (^2 _ ^2)^(1 _ + ^2^(1 _x)-0 + ief 

(A.3) 

^ _dH 2(1 - x)x{p'^ - I?) 

(27r)4 [p + (p2 _ ^2)3,(1 _ + ^2^,(1 _x)-0 + ief 





where O = ax — P{1 — x). Integrating over the loop momentum and making substitution 

y ^ u: 

y x{l-x)y 



-20- 



we obtain the DR's 



1 oo 



dx 



duo 



o 



1 — X — /i^ 

Lo — fi"^ {p^ — Lo + ie) 



x{l — x) 

oo 



= C'fo j dx j duj 

o 

x{l~x) 

Using the fohowing definition of X functions 

1 

X{u!;a,P) = J dxQ ( u> 





uj — pfl {p^ — uj + ie) 



O 



x{l — x) 
O 



r ( 

X\ (lu; a, P) = / dxxQ to — 

J V ^(1 





we can rewrite DR( A.5 ) into the more famihar form: 

.2 r 



(a; — ^J?){p'^ — to + ie) 
X{u;;a,P){p^ -fi^) 



{ui — iJ?){p'^ — oj + ie) ' 
where we made use of the relation 

a, /3) — a, /3) = /5, a) 

Exphcit integrations in ( |A.6D give the fohowing formulas for X functions 



Xi (w; a, (3) 



A^/^(a,cu,/3) 
2uj 



1 + 



13 — a 



UJ 



(A.5) 



(A.6) 



(A.7) 



(Ai 



(A.9) 



where A is the triangle Khallen function A = (w — a + — 4u;/3. 

After the explicit introduction of the prefactors Ca,b we can write down the absorptive 
part of self-energy functions a, h: 



-9 



rfr^Xi{uj\m^ ,m ) + mr^ j daa^{a)Xi{oj;m^,a) 



+ 



(47r)2 
oo 

Tf J dPp^{P)X,{u;;P 



m 



+ J dp J dap^{P)av{a)Xi{uj; P,a) 

Am? (m+ra^)'^ 



(A.IO) 
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(47r)2 



daas{a)X[u); m^, a) 



Am? 

oo oo 

+ j df3 j dap^{l3)as{a)X{uj;a,l3) 

im? (m+m^)^ 



(A.ll) 



As follows from the properties of function X the various terms in ( A.10),( A.11| ) start to be 
nonzero from different values of uo. This knowledge is particularly useful when calculated 
numerically. Indeed after the inspection of X{uj] a, (5) we can find the threshold at a; = 
V^). For instance the subtresholds values (m^ + m)^; {2m^ + m)^; 9m^; (m^ + 3m)^ 
correspond to the terms in Eq. ( A.ll| ) at given order. 



A. 2 Yukawa pseudoscalar self-energy 

Let us find the DR for self-energy 11. The loop integral has the form: 

dH ^ ^5[^4a)+cT,](a)75[(F^)^T.(/3) + ^s(/?)] 



(2vr)4 
dH 



Tr- 



- a + ie){{l - pY - 13 + ie) 



{{l-py^P + ie) 



1 



+ / dx 





(-g - x))q^^,(a)^7t,(/3) + as{a)crsiP) 

P + p'^x{l — x) — O 



(A.12) 



where the shift I I +p{l — x) is made at the second integral and O = ax — /?(! — x). For 
purpose of brevity we omit the spectral integration over the variables a and /?. 



The renormalization (2.7) proceeds by direct subtracting of the first two term in Tay- 
lor expansion of unrenormalized self-energy ( A.12| ). To make this explicitly we introduce 
shorthand notation U{w,l,0) = — w^{l — x)x — O + ie) . After a little algebra the 
renormalized quantity corresponding to ( A.12| ) can be evaluated as 



n(/i;/ 



4i 



dH 



dxdy 



+ 



2x(l 



UHfi,l,0) 



U^{p,l,0) U^{^^,1,0) 

- x) fi'^{l - x) ^ 2x(l - x)^fi^{p'^ - fi^) 



u^p,i,o) u^p,i,o) ' uHp,i,o) 

Matching terms in the first line of ( [A. 131 ) together we can rewrite them into the familiar 



DR 



(4vr) 



duj 



(w — p?Y{p'^ — uj + ie) 



2\2 



(A.14) 
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where function X was already introduced in ( A.9 ) 

The terms in the second brackets of the second hne of ( A.13| ) can be matched together 
in the following fashion: 



± 

/ 



dy 



2(r - - x) 



{Uip,l,0)y + Uifi,l,0){l-y))^ 

2, ,2/ 2 ''-^ ' 1 



+ 



(2x(l-x)V(p'-/x')) 



UHfi, I, O) iU{p, 1, 0)y + U{fi, 1, 0)(1 - 



dy 



+ 



7, I T)2 ,,2 _|_ /£! _ O 

y\P - P' + y x{l-x)y 



{p' - 



O 



II Z ^ \ I nZ /i2 _l_ A* w 

^ ~ x{l-x)y [P ~ ^ + Y ~ x(l-x) 



O 



(A.15) 



Adding omitted prefactors, integrating over loop momentum and making the substitution 
( |A.4| ) in ( |A.15| ) gives rise to the DR: 



(47r)2 



2\2 



duj 



(X(a;;a,/?)-Xi(cu;a,/3)Mp^ -/x^) 
{u — — w + ze) 



(A.16) 



Putting all together and using the relation between X and Xi we can finally write down 
the DR for self-energy 11: 



n(/i;/) 



doj 



2\2 



(u;^ — — u) + ie) 

da / dp [{—aay{a)av{f3) + as{a)as{f3)) X{u!;a, (3) 



(47r)2 



(A.17) 



For the purpose of numerical solution it is necessary to extract the singular parts of 
SG's explicitly. After some trivial manipulation the appropriate formula for reads 

2 



a_ 

27r 



1 



Am? 



to 



+ 2rjm 



oo 



+ 



(ni+m^)2 



da 



dp { -{to - a- P)aTj{a)ay{P) 



+ as{a)(Js{P)) X{LO]a, 



(A.18) 



noting that the first line of ([A.18| ) corresponds with the one loop perturbative contri- 
bution. 
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B. Quenched, rainbow, Landau gauge QED 



It is shown, that with a massless photon the fermion self-energy function b satisfies DR 
(2.18) while the function a is exactly zero. In the end of this section we review the results 
for the case of massive photon. These are used in the numerical study described in the 
main text. 

First we will deal with massless photon case (A = 0), where our approximate electron 
self-energy 



requires one subtraction due to the presence of logarithmic ultraviolet divergence. As it is 
sometimes usual, we fist regularize (B.l), then the subtraction will follow. For this purpose 
we use the old-fashion Pauli-Villars regularization technique which leads to the following 
regularized result: 

^a(^) = ~J^^ j j ^^t*^^'''^") ^(4x - 2) + 3a-s(a)) , 



, , -p'^x(l — x) + ax + A'^il — x) 

X In ( -^-^ ^ -— ^ + ie 

— p^x(l — X) + ax 



^ 2ay{a) ^p^x^(l — x) 



(B.2) 



p'^x{l — x) + ax 

where A represents Pauli-Villars regulator. Making the aforementioned subtraction gives 
the renormalized functions a, b 



a{fi;p 



da 



ay{a)e'^ 



(4vr) 



+ 



2p2 



dx{4:X - 2) In 



2/i2 



-/U^(l — x) + a 
-p2(i — x) + a 



+ ie 



-p^(l — x) + a —^^(1 — x) + a 
1 



6(^;p^ 



da 



2>as{a)e 
(4^)2 



dx In 



-11 [l — x) + a 
-p^{l — x) + a 



+ ^e 



(B.3) 



where /i^ is space-like renormalization scale. The per-partes integration of logarithm terms 
gives 



6(/x;/ 



da 



3(Ts(a)e^ 
(47r)2 



dx 



— ^ I da I duj 



(1-x) {p2-^ + ^e) 
3as{a){l - ^){p' - fi') 



{uj — ^'^){p'^ — uj + ie) 



(B.4) 
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Thus absorptive part of self-energy -Kp is given by : 



rm [ 1 ) + / daas{a) ( 1 



(B.5) 



Note here, that the Pauh-Villars regularization technique was chosen for convenience only. 
It is not so difficult to use, for instance, the dimensional regularization technique. Of 
course, making direct algebraic subtraction of unregularized a, b is also possible. All these 
purely technically different approaches leads to the same results, that is the subject what 
we exactly understand under the statement " regularization independence" . 

When the photon propagator Gq is changed by the introduction of small mass param- 
eter A 



(B.6) 



A;2 - A2 + ie 

then the previous results are slightly modified. Repeating the steps as above it leads to the 
same form of dispersion relation ( ]2.18 ) but with different ps,Pv j both of them are nonzero 
now. We get for them 



Ps{^) = ~3 (^) |mX(u;, m^, A^) + j (iacrs(a)X(a', a, A^) | 
p^{uj) = W X(u;,m2,A2) - 2Xi(cj,m2,A2) 



A2 

Xi{uj,a,\^) - 2Xi{uj,a,Q) 
A2 



+ [u) — m ) h / daas[a) X(uj,a,A ^ 



2Xi{ijj,a,\^) + {uj -a) 



(B.7) 



Stressed here, that the limit A — > can be safely performed an leads to the result (|BT 
and p^ = 0. 



C. Numerical details 

The unitary equations have been solved by the method of numerical iteration. To achieve a 
reasonable accuracy we should carefully perform the principal value integrations labeled by 
[cj, p], ([c7, p\) in the equation for fermion (pseudoscalar) weights. For instance the numerical 
P. integration 

'y{s)p{x)'-^ + a{x)p{s) 



[a * p] = p. / dx- 



s — X 



(C.l) 



thresh. 



is proceed by the following way which is based on the exact relation = / dx ^ . Hence 



we can write for P. integral in ( |C.l| ): 



oo oo 

p. I dx^ = P f ^ - 



a — X 

c c 



a2 — x2 



'a + c 



(C.2) 
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where 

f{x) = cr{a)p{x) 2 + cr{x)p{a). 

The right hand side of the identity (|C.2D is particularly useful when evaluated numerically 
i.e., when J ^ S. 

Using some contemporary PC machine the criterion o"^ ~ 10~^^ can be achieved, 
noting that the typical CPU time is about 10^ s for a grid number of several hundred mesh 
points. Here, o"^ represents the error between the solution of n times and n — 1 times 
iterated unitary equation i.e., 

2 _ liPn ~ Pn-l) o\ 
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Figure 1: The spectral functions for Yukawa model calculated for the values 0.1; 0.2; 0.3 of coupling 

strength A = The lines with the threshold (m + m^)^ = l.lb'^m? correspond with the Dirac 
(negative lines) and scalar (positive lines). The Lehmann weights for pseudoscalar propagator have 
the appropriate thresholds at Arn? and they are always positive. 
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Figure 2: The self-energy of pseudoscalar meson Q(p^) = m'^{fi) + ni?(/i^,p^) for various coupling 
strengths A of Yukawa interaction. The only real part is displayed above the threshold 4rn^. 
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Figure 3: Yukawa fermion functions A{p'^), B{p'^) at space-like regime of momentum. 
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Figure 4: Yukawa fermion fimctions j4(p^), B{p^) at time-like regime of momentum. 
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Figure 5: The renormalization scale invariant dynamical mass of Yukawa fermion. The renor- 
malized mass to(0) is fixed at /i = such that the pole mass is m/ = f for all the couplings 
A = 0.1;0.2;0.3 
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Figure 6: The dynamical mass of the electron in massive photon QED as they have been obtained 
by solving modified unitary equtaions. The various lines are labeled by the mass of the photon A 
which is written in the units of an on-shell electron mass. All solutions have the coupling value 
a = 0.6, renormalization point /i^ = —10^, and renormalized mass is = 400. The results for 

two smallest A are not distinguishable. Thin long dashed line labeled by the letter E represents the 
solution obtained in Euclidean formalism. In this case, photon was exactly massless. 
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Figure 7: The electron dynamical mass M{p^) for space-like regime of momentum. Six dotted lines 
represent the Euclidean solutions of momentum DSE. They are labeled by the value of coupling 
constant at given order (from up to down). The solutions of unitary equations (dot dashed lines) and 
the ones of modified unitary equations (solid thin lines) are added for the coupling a ~ 0.2; 0.4; 0.6. 
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Figure 8: The electron dynamical mass M{p'^) = B{p^) for time-like regime of momentum and 
renormalization choice m(0) = 1. The thick lines represent Euclidean solutions obtained by principal 
value integration of the Fukuda-Kugo equation. The thin lines correspond with the solutions of 
modified unitary equations (labeled by UE). 
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Figure 9: The same as in tlie previous figure, but for stronger couplings. The results are obtained 
from solution of Fukuda-Kugo equation only. Signal for confinement is apparent for the coupling 
a > 0.915. The lines represent a good estimate of true FKE solution for subcritical couplings but 
they must be rather exact for the supercritical couplings solutions (for the meaning of e see ([4.3|)). 
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Figure 10: Space-like solutions for electron propagator. All solutions have the renormalization 
point fi^ — —100, and renormalized mass m(/i) = 10 and couplings is varied from 0.1 to 8.0. Note 
that there are several zero-crossing when the coupling is large enough. The solutions of modified 
unitary equations (dashed thin lines which are closed to their Euclidean counter-partners) are added 
for the lowest couphngs a = 0.1, 0.2, 0.3, 0.5. 
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Figure 12: The absolute values of absorptive part of electron self-energy. In order to obtain the 
result for a = 1, we used the infrared cut-off value 0.2m in the modified unitary equations treatment 
(m is a pole mass, here always m ~ m(— 100) = 10, for instance m = 14 for a = 0.8. 
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